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Abstract 



Including nucleon-nucleon correlations due to both Fermi statistics and 
nuclear forces, we have developed a general formalism for calculating the 
neutral-current neutrino-nucleon scattering rates in nuclear matter. We de- 
rive corrections to the dynamic structure factors due to both density and spin 
correlations and find that neutrino-nucleon scattering rates are suppressed 
by large factors around and above nuclear density. Hence, in particular for 
the and u T neutrinos, but also for the v e neutrinos, supernova cores are 
more "transparent" than previously thought. The many-body corrections in- 
crease with density, decrease with temperature, and are roughly independent 
of incident neutrino energy. 

In addition, we find that the spectrum of energy transfers in neutrino 
scattering is considerably broadened by the interactions in the medium. An 
identifiable component of this broadening comes from the absorption and 
emission of quanta of collective modes akin to the Gamow-Teller and Giant 
Dipole resonances in nuclei (zero-sound; spin sound), with Cerenkov kinemat- 
ics. 

Under the assumption that both the charged-current and the neutral- 
current cross sections are decreased by many-body effects, we calculate a set 
of ad hoc protoneutron star cooling models to gauge the potential importance 
of the new opacities to the supernova itself. While the early luminosities are 
not altered, the luminosities after many hundreds of milliseconds to seconds 
can be increased by factors that range from 10% to 100%. Such enhancements 
may have a bearing on the efficacy of the neutrino-driven supernova mecha- 
nism, the delay to explosion, the energy of the explosion, and the strength and 
relative role of convective overturn at late times. However, the actual conse- 



quences, if any, of these new neutrino opacities remain to be determined. 
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I. INTRODUCTION 



In core-collapse supernovae, neutrinos are arguably the engines of explosion and the 
best direct probe of its internal workings To understand the supernova phenomenon, 

a theorist requires knowledge of the equation of state of nuclear matter, stellar evolution, 
general relativity, statistical physics, and the techniques of radiative transfer and hydro- 
dynamics. However, the opacities and sources of neutrinos of all the six known neutrino 
species take center stage in this context in which the neutrino "optical" depth down to the 
protoneutron star core varies from 10 3 to 10 6 , neutrino luminosities can approach 10 54 ergs 
s _1 , and photons, the traditional agents of radiative energy transfer in astrophysics, are 
profoundly trapped. This seems to be a unique role for neutrinos in the universe. Since neu- 
trinos are abundantly produced, yet are the most mobile of a collapsed core's constituents, 
their microphysics determines the outcome of core collapse, neutron star and black hole 
formation, nucleosynthesis at and beyond the iron peak pH|jrT[ |, early pulsar kinematics |T2| , 
and, perhaps, the asymmetries observed in supernova debris clouds |13| . 

In the neutrino-driven mechanism, neutrinos liberated from the core heat matter in the 
inner stellar envelope (the outer protoneutron star envelope) as they emerge. If the total 
heating rate is sufficient, the envelope becomes unstable and is ejected. It has recently been 
shown P|-P|,P^| that neutrino-driven convection in this inner mantle can make this process 



more efficient (but see |15]). Hence, neutrinos mediate the energy transfer from the core to 
the mantle and their luminosities and spectra are crucial to the viablility and character of 
the explosion mechanism. These in turn are a function of the neutrino-matter opacities. 
While at low mass densities (p < 10 12 gm cm -3 ) neutrino-matter cross sections are well- 



understood and characterized |16|-|T^,^] , at the higher densities achieved in core collapse 
the inter-particle spacings are smaller than the wavelength of the ambient neutrinos and 
many-body effects and particle-particle correlations must be taken into account. Ion-ion 
correlations due to Coulomb interactions have been addressed in the supernova context (cf. 

and their effects, though interesting, have been shown to be small and transient |20| . 
However, near and above nuclear densities (~ 2.6 x 10 14 gm cm -3 ) neutrinos with energies 
below ~1000 MeV "see" many nucleons and collective effects must be considered. In fact, 
these effects can radically alter the opacity of nuclear matter to neutrinos, since the nuclear 
force introduces strong correlations near and above nuclear densities. 

Many-body effects have been discussed before in the protoneutron star and supernova 
context (cf., [pTH 25|1 ) , but the results were either too approximate or of too-limited scope. 
Pauli blocking by final-state nucleons, a type of correlation, has been a component of super- 
nova and protoneutron star thinking for some time P J26"1 , |2"7|1 , but the approximate treatments 
employed at high densities can be off by large factors (§VI.B). Furthermore, in the case of 
neutrino-nucleon scattering, the assumption that the scattering is elastic |28[| and that the 
nucleons are stationary can be shown to result in errors in the energy transfer of as much as 
a factor of ten. A consequence of this is that energy equilibration of and v T neutrinos is 
brought about not by the v + v — > e + + e~ process or neutrino-electron scattering, but by 
neutrino-nucleon scattering. This result, one by-product of this study, has also been found 
by Hannestad Sz Raffelt ]29| and in the pioneering work of Reddy et al. ||30| . The latter 
have conducted the most complete exploration of high-density neutrino opacities to date, 
taking care to calculate them in a manner consistent with a reasonable nuclear equation of 
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state (see also |pT| , |32| ). Some of our results and formalism recapitulate theirs. However, 



Redely et al. |30[ neglected particle-particle and spin-spin correlations and focussed on the 
kinematic effects of nuclear interactions. We find that the effect of many-body correlations 
on neutrino-matter cross sections can be quite large. Nuclear matter is more "transparent" 
to neutrinos than we had heretofore imagined. 

Neutrinos interact with nucleons via both charged-current and neutral-current processes. 
The charged-current processes predominate for the electron-type neutrinos by a few factors, 
while the neutral current predominates for and v T neutrinos. The latter carry away 
more than 50% of the binding energy of the neutron star. We have found that at the 
high densities achieved in collapse, both the charged-current and the neutral-current cross 
sections are dramatically affected by nuclear correlation effects. In this paper, we focus on 
neutral-current neutrino-nucleon scattering, both the Fermi and the Gamow-Teller parts, 
and defer a discussion of the charged-current many-body effects to a later paper. Our 
formalism is good for any degree of nucleon degeneracy and fully incorporates the effects of 
reaction kinetics, Pauli blocking, and correlations due to interactions. Energy transfers both 
to and from the nucleons are consistently included. One product of this paper is a formalism 
for calculating in full the dynamic structure functions for neutrino-nucleon processes. 

The effect of increased transparency at high densities on the neutrino-driven mechanism 
of core-collapse supernovae is not yet clear. Supernova theory is notorious for its mitigating 
feedbacks and false leads. The early (,$,500 milliseconds) neutrino luminosities and spectra 
depend upon the material in the outer shocked mantle of the protoneutron star, where we 
think we understand neutrino cross sections. As a consequence, we don't expect that our 
understanding of this early phase will be altered by the new opacities. Even if it could be 
shown that neutrino cross sections at the lower densities around the neutrinospheres and 
near the shock were in some non-trivial sense different, increased transparency would lead 
to an increased rate of collapse of the mantle, which, in turn, would lead to higher densities, 
which would partially quench the effect. The corrections to the scattering rates we address in 
this paper are more relevant at later times (> 500 — 2000 milliseconds), for it is then that the 
emergent luminosities are powered by the energy in the dense core. At this time, since most 
of the core pressure is derived from the cold, stiff nuclear component, and not an ideal gas of 
non-interacting nucleons, the effective specific heat is not negative, but positive, and cooling 
does not lead to much of a density change. Therefore, the effect of a decrease in opacity will 
not be partially cancelled by a density increase and the emergent luminosities will reflect the 
full suppression effects. Since the neutrino luminosity is the agency of explosion, this could 
be important for a model in which the delay to explosion is not short. There is an important 
caveat: accretion of the infalling envelope is a major power source and its continuance can 
mask the effects of the alterations we find in the neutrino opacities at depth. Accretion 
will be less important for the lower-mass progenitor stars (8 — 13 M Q ), which have more 
tenuous envelopes. In addition, the neutrino signals observed from SN1987A |],[/]] constrain 
the magnitude of any alteration in our standard model, as does the fact that there is a fixed 
amount of energy to be radiated from any given protoneutron star. In this paper, to gauge 
these effects and to begin the debate on these issues, we perform an idealized protoneutron 
star cooling calculation. 

In §11, we explain the physics of the correlation effects and provide a simple, single- 
channel, model that makes clear the character of the results. In §111, the general formulation 
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is presented, as is an aside on the low-wavelength limit. In §IV, the complete correlation and 
dynamic structure function formalism is derived and discussed and a straightforward nucleon 
potential model is described. In reality, the solution of the problem of the nuclear equation 
of state goes hand-in-hand with the calculation of many-body correlations. Since we are not 
proposing here to resolve the former, we are satisfied with a model that captures the essence 
of the nuclear interaction. Different nuclear equations of state will yield quantitatively, 
but not qualitatively, different results. In §V, we outline our multi-channel formalism for 
calculating //-nucleon scattering structure functions and in §VI we present the results of our 
calculations for various temperatures, densities, and neutrino energies. In §11 through §VI, 
we start from the general expressions of statistical mechanics that encompass the effects of 
the medium on the neutral-current neutrino reactions and proceed to calculations of the 
effects. In §VII, we present the results of a suggestive series of calculations that might point 
to the astrophysical import of the new effects and in §VIII we summarize our conclusions. 



II. THE ESSENTIAL PHYSICS 

We can introduce, and qualitatively explain, much of what is to come in a simplified 
system in which we have just one species of non-relativistic nucleon with mass m, which is 
coupled to neutrinos through only a vector neutral current. We take the weak-interaction 
Lagrangean density to be, 

G 

C w = -j=ij) v (x)(\ - lb)lo^u{x)n{x), (1) 

where n(x) is the density operator for the nucleons, Gw is the weak coupling constant, and 
the other symbols have their standard field-theoretic meanings. In the following, we employ 
a system of units in which h — c — ks — 1. 

We write the differential rate of neutrino scattering p x — > p 2 in the medium as 

<]2{ = (AnY 1 G 2 w (E 1 -uf[l-f u (E 1 -co)]A 00 (q,co)S(q,u J ), (2) 



dtu d cos 9 



where f u is the occupation function for neutrinos, E\ is the energy of the incident neutrino, 
A is the neutrino trace, 

A"" = ^E x E 2 )- x Tr[i x (1 - 75)7% 7"(1 - 7s)], (3) 

uj is the energy transfer (|pi| — IP2I) to the medium, and q is the momentum transfer 
(|Pi — P2I), related to u> and pi through the neutrino scattering angle, 9, by 

q = [Pi + (Pi - - 2pi(pi - u) cos#] 1/2 . (4) 

For the case of free nucleons, the nucleon part, S(q,u) (the dynamic structure function) 
is given by 

S(q,v) = 2f ^/(|p|)(l-/(|p + q|))27r% + 6 p -6 p+q ), (5) 



5 



where /(|p|) is the Fermi-Dirac distribution function and e p is the nucleon energy. Eq. (0) 
is Fermi's golden rule with final-state nucleon blocking. The energy delta function expresses 
energy conservation and momentum conservation, comes from the matrix elements of the 
nucleon density operators < p'\ J ci 3 xe jq ' x n(x)|p >, and has already been integrated out. 
From the 5 function in (|5]) we see that a typical energy transfer u from the medium to the 
neutrino is of the order of q times the thermal velocity, (T/m) 1//2 . Thus, in the limit of heavy 
nucleons, when we integrate the differential rate (0) over a range of ujs the other factors in 
the integrand can be evaluated at uj = 0. We can express this limit as, 

{27i)- 1 S{q,u) -> {2tt)- 1 5{uj) f dJ S{q,uj') = 5{uj)S{q), (6) 



where S(q) is the static structure factor. In §VI, we show that at the high densities and 
temperatures achieved in the supernova context the uj = (elastic) limit is not particularly 
accurate (see also fl30"f ). 

When we turn on interactions among the nucleons, (|J) is replaced by, 

S(q,oj) = 27TZ- 1 J2 e _/3 ^ / rfW qx < j\n{x)\k >< k\n(0)\j > Sfa + Ej-Ek), (7) 

where j and k are energy eigenstates, now of the whole medium, Z is the partition function, 
T is the temperature, and /3 = ~. If the medium is comprised of heavy enough nucleons, we 
need only the static structure function defined in (^), then given by 

S(q) = Z- l Y. e ~^ E] J d 3 xe i(l - x < j|n(x)n(0)|j > . (8) 

j 

In other words, S(q) is merely the Fourier transform of the thermally-averaged density- 
density correlation function. This is the classic result that scattering off of a medium is 
in reality scattering off of the fluctuations in that medium. Also of interest is the long- 
wavelength limit, q — > 0, justified when the neutrino wavelength is much bigger than the in- 
terparticle separation.^] Statistical mechanics provides two useful and equivalent expressions 
for the long-wavelength limit, 5(0), the first |33j in terms of the isothermal compressibility 
of the medium Kt ( : 



dlogV I 
OP \T) 



S(0) = n'r'Kr = f&, (9) 







where K Q is the ideal gas compressibility and n is the average nucleon density, and the 
second {ibid., their eq. 114.14) in terms of the derivative of the density with respect to the 
chemical potential of the nucleons, /i, 



x We emphasize that the limits discussed in this section are presented both for pedagogical reasons 
and because they provide a boundary condition for the later work in this paper, in which we require 
neither limit. Under the actual conditions that prevail in the supernova core the conditions for the 
limits to be applicable are marginally satisfied, at best. 
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In the ideal gas limit of no correlation between particles, eqs. (|9]) and ( |T0D show that S(q) is 
simply equal to the number density, n, as expected from eq. (|8|) and eq. (|5|), without blocking. 
Eq. (|J) reveals that if Kt is small because the matter is stiff, in the long-wavelength limit 
the neutrino-matter cross sections are suppressed. When we replace eq. ([I]) by the complete 
standard model form, including the axial-vector current and nucleon isospin, we shall require 
separate correlation functions for the neutron and the proton, as well as for spin correlations. 
These depend upon susceptibilities that are different from the compressibility, but we shall 
find suppression in these terms as well. 

Eq. (|lOl) , equivalent to eq. (§J) by a thermodynamic identity, is a powerful result of great 
generality. In standard approximation schemes for the many-body problem, the distribution 
function for a nucleon species is given by the Fermi-Dirac distribution in which the chemical 
potential [i is replaced by fx — v(n), where v (n) is the average energy of interaction of the 
nucleon with the other nucleons and is a function of the density. Thus, the density is given 
implicitly by 

n = 2 J Jl!L[1 + e /nb7(2™)-/^)]]-i. (n) 

The expression ( |TT| ) holds in the Hartree approximation; it holds in approaches that 
introduce mean meson fields instead of nuclear potentials; it holds in the Landau Fermi 
liquid theory (FLT), subject to the proviso that we use only results in which the derivative 
of the potential v (with respect to the n) enters; and it holds in approaches using the Skyrme 
potential. 

Differentiating: (|IT|), we can solve for g and S(0), 

S(0)=(3- 1 ^- = h^)[l + h(^}-\ (12) 

where 

f d 3 p e P[p 2 /(2m)-Li+v] „ 

h{fl) =2 J ( 27 r)3[l + e^V( 2m )-^]]2 = 2 J J2^f h ^ 1 ~ (13) 

and fi{p) is the Fermi-Dirac function, but with the chemical potential displaced by v. If we 
regard particle densities as inputs to our calculations, then the displacement of the chemical 
potential by the nuclear potential is irrelevant, since the same difference, fi — v, enters 
the calculation of the density in terms of the chemical potential. Thus, the numerator of 
(|l2l) contains no more than the familiar Pauli blocking effects (for the case q = 0); the 
denominator contains all of the effect of the interactions. 

As an example, consider a two-nucleon potential V(r). In the Hartree approximation, 
the average potential seen by a single nucleon is given by v = nil, where U = J d 3 xV(x), 
and fll2|) becomes 

S(0) = h(fi)[l + h(fi)Ur\ (14) 
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the potential providing an enhancement, if negative, and a suppression, if positive. The 
latter is the case for high-density nuclear matter. 

Extensions of the above considerations to include isotopic spin and the axial-vector weak 



hadronic current were given in reference ||34 | . Closely related considerations were given in 



reference [2TJ for the case of degenerate neutron matter. Both references [M] and |2T 
concluded that nuclear interactions cause a big decrease in the Gamow-Teller part of the 
neutrino opacity in the regions considered in the respective works. However, in addition 
to being fragmentary in the domains that were covered, these works did not address two 
important issues: 1) the errors in doing the integral over u) by assuming that the neutrino 
parts of the matrix element and phase space factors are independent of u over the region 
of dominant contribution and, 2) the applicability of the q = limit. These limits are the 
least justified for the case of trapped electron neutrinos in the early dense core, where the 
momentum transfers in neutrino scattering are comparable to the inverse particle spacing. 

In the present work we address these deficiencies by using methods that incorporate 
all of the physics of the above discussion and previously cited papers, but which calculate 
the dynamic structure factor as a function of (q, u) and do the integrals over the neutrino 
variables, without further approximation. We also aim to better systematize the problem 
by keeping the connections to textbook many-body formalism Q as clear as possible. 

III. GENERAL FORMULATION 

To address in complete form the reactions that a neutrino undergoes in a medium we need 
to calculate the appropriate thermal averages of the weak-current operators for the particles 
that comprise the medium. We define as the weak neutral current of the particles in the 
medium, other than the neutrinos, and we assume conditions under which all species, except 
the neutrinos, are in statistical equilibrium. The properties of the medium are embodied in 
the function W^ v 

Wp,(q,u) = F[j^3v\q^ , (15) 
where we define the correlation functional T q ^ of two Heisenberg operators 0\ and O2 as 

F[0 1 ,0 2 ] q ^ = Z- 1 J d i xe- l ^e luJt Tr[e- l3{H - E ^ ) O 1 (-K,t)O 2 (0,0)}, (16) 

where Z is the partition function. We can write the differential rate for neutrino scattering 
in terms of W^ u (q, u), 

, f F a = (4tt 2 )- 1 G'^(S 1 - uf[l - UEr-u^iq^W^)- (17) 
dud cost) 

Note that the leptonic spinor-trace and phase-space factors are as in the golden rule 
formula (§). The distribution functions for initial and final nucleon states, and the golden- 
rule delta function, have all been subsumed in W^q.u). When the medium is interacting, 



2 In what follows, we will use a number of formulae from Fetter and Walecka hereafter referred 
to as FW. 
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(|17|) can be taken as the definition of the correct answer, as it is in the form that ab initio 
thermal field theory gives for neutrino transition rates. Alternatively, when a complete set 
of intermediate energy and momentum eigenstates are inserted between the two operators 
in fll6|), the d A x integrals give the delta functions, the two matrix elements give the square 
of the T matrix, and the thermal factor weights them appropriately, as in ([?]). 

We take the medium to be composed of protons, neutrons, electrons and neutrinos and 
we deal only with the neutrino scattering from the nucleons in the medium.^ 

For the nuclear contribution, we assume that the nucleons remain non-relativistic to 
write the space components of the function W^ u in the form, 

Wi, = W 0>i = Wij(q, w) = W A {q, oj)8 itj + W B {q, u)q iqj . (18) 

We define Wy = Wo,o- The lepton trace can be expressed in terms of these functions as 
follows, 

W^tsT = (l+cos0)Wy + (3-cos0)M^ + w 2 (l+cos0)WB, (19) 

where, as before, 9 is the scattering angle, Wy represents the Fermi (vector) term, and Wa 
is the Gamow-Teller (axial-vector) term. 

The Wb term in eq. fliHf) vanishes in the absence of interactions among the nucleons. 
In our later parameterization of the forces, it can become different from zero only when a 
tensor force is added to the conventional four central forces. We have directly evaluated the 
contribution to the opacity from the tensor force coming from single neutral-pion exchange 
and find, for degenerate matter around nuclear density, Wb ~ 2 x 10 _5 W / y i(MeV) _2 . The 
tensor terms give an additional opacity that is a small fraction of that due to the axial-vector 
terms (Wa), and we drop them in what follows. 

The non-relativistic limits of the standard model neutral current to be used in (|T5|) can 
be expressed in terms of the neutron and proton density operators, n n and n p , and the z 
component of the spin-density operators, n^ p = V'n.p^V'n,?) 

jo (a?) = Cy ] n p (x) + Cy ] n n (x) 

j 3 (x)=g A (n p 3 \x)-nW(x)), (20) 

where eft = 1/2 - 2 sin 2 9 W , C { y ] = 1/2, g A = -1.26/2, and sin 2 6 W = 0.23. By rotational 
invariance, we need only the correlation functions of the z component of the spin, but, since 
the medium is not invariant under isospin rotations, we use separate neutron and proton 
currents. 

For the vector (Fermi) part, we now obtain, 

W v (q,u>) = (C^) 2 S pp (q,uj) + 2C^C^S pn (q,uj) + (C^) 2 S nn (q, u), (21) 



Of course, neutrino-electron scattering through the neutral-current couplings must be taken 
into account in the determination of the total opacity, and for the case of electron neutrinos the 
charged-current interactions with nucleons and with electrons must be used as well. 



9 



where the structure functions are defined as S pp = J-[n p , n p ], S pn = J-[n p , n n ], and S nn = 
F[n n ,n n \. 

For the axial-vector (Gamow-Teller) part, we have 

W A (q, u) = g 2 A (S^(q, u>) + St(q, u>) - 2S*(q, u)) = g 2 A S A (q, u), (22) 

where S£ = F[nf\nf\ S^ n = F[nf , n^} , and S* n = F[n£\nV>]. The purpose of the 
remainder of this paper is to calculate the six structure functions. 



A. An Aside on the Long Wavelength Limit 

In the long-wavelength limit, we can calculate the six structure functions quite sim- 
ply and intuitively. We introduce the average densities, n*, for the four species, fti = 
(Vol) -1 J d 3 xrii(x), where the index % runs over the values p* , p^-, n* , n^, and we form 
the combinations, 

n Py n = Tlp.n "I" ^p,n 

^p,n = ™p,n ~ ^p,n- (23) 

We also introduce separate chemical potentials for the up and down spin states for the 
two species and introduce the notations, 

i T I ^ 

^P,n = Q \^p,n ' ^p,n) 



Then, as shown in ( |T0"D and reference |34 l, we can express the q — > limits of the static 
structure functions as follows: 

Spp(0) = (3 ; S nn (0) = P' 1 ^ 

s pn (o) = s np (o) = 

^(o) = r 1 || ; st(o) = p^ 9 4 

Spi(0) = S^ p (0) = (3-^. (25) 

As in the one-channel case treated in §11, we take the distribution function for a nucleon 
species, i, to be given by the Fermi-Dirac distribution in which the chemical potential, /ij, 
has been replaced by //j — u,, where V{ is the average energy of interaction of the nucleon 
with all other particles. Thus, the density is given by 



71: 



d 3 P r-i , D /3(p 2 /(2m,)-p l +v t ) ] -l 



[1 _|_ gPtP /(2m;-Ati+tH;j-i_ ^6) 



In general, the potential, Vi for species % is a function of the densities for all four of the 
species, fij. If we know the functional dependence of the v's on the UjS, we can solve for 
the long-wavelength limit of the Fermi static structure functions, S, by differentiating the 
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four equations (|26|), with respect to the four chemical potentials that we have introduced. 
We give the solution for the two combinations of structure functions that enter our rate 
calculations. In (|21|), with sin 2 6w = 1/4, we have only the neutron contribution, 



where 



Srt 



d v = 1 + ^-h(n n ) 
on n 



P' 1 — — = h(fx n )/d v , 



h(nn)h(fi p ) 



dv p hi 



(27) 



(28) 



In (| 

dvi I drij 

by the potentials we use. We then obtain, 



, integrated over u, we obtain a simple form, if we choose potentials such that 
dvj/dfii and such that v p <\ — v p ± + u n j — v n ± = 0, two conditions that are fulfilled 



S A (q) 



(q) 



h(fi p ) + h(fi n ) 
l + v a (h(n P ) + h(fx n ))' 



(29) 



where Ua = -[^ - g^][v p1 - v pl - v n1 + v nl }. 

We present these static, long-wavelength results (eqs. |27| and |29|), which depend on ef- 
fective single-particle potentials, because they show a strong similarity in form to the more 
complete equations that will be developed in the next section and because, in most regions, 
they give numerical results that are not greatly different. However, it should be pointed 
out that any program that begins with potentials that fit nuclear data and calculates the 
ground-state properties of nuclear matter (c/., [0,0) is capable, if subjected to the right 
sets of constraints, of directly determining the low-temperature values of the static structure 
functions enumerated in fl25|), without recourse to the assumptions that led to eqs. ( p7|) and 
(p5p. Implementation requires use of the multichannel analogues of the connection between 
the one-species structure function and the bulk modulus (|9]). The required numerical exper- 
iment involves constraining the system to have different expectation values of proton spin 
and neutron spin densities (in the z direction), calculating the constrained ground-state 
energy as a function of these densities (as well as the particle densities), and taking combi- 
nations of second derivatives of this energy with respect to the densities, evaluated finally 
for the values of the densities in the true ground state. We strongly recommend that groups 
that do nuclear matter calculations carry out these steps. 



IV. COMPLETE CORRELATION FUNCTIONS 

In the field-theoretic formalism for the quantum mechanical many-body problem there 
are well-formulated perturbative techniques for directly calculating correlation functions, 
without restricting to the equal time (u = 0) and q — ► limits. In Appendix B, we show 
that the limiting results sketched out in the previous section by calculating the derivatives, 
1^-, are exactly what we obtain, in the same limits, from the sum of the "ring" approximation 
graphs in this many-body formalism. Therefore, we must sum at least the ring graphs to 
recapture the long-wavelength limit of the static structure function. In so doing, we obtain 
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the full dynamic structure function S(q, uj) at the same level at which we find the limiting 
forms in §111. A. 

To follow a graph summing approach we must replace the "correlation functional" of two 
(bosonic) operators defined in (|T6|) , JF, by a retarded commutator form, T, 

j r i0 1 ,0 2 ] q ^ = ~ J d 4 xe-^e^ t Tr[e-^ H - s ^[O 1 ( X ,t),O 2 (0,0)}]9(t), (30) 

where here 9(t) is the Theta function. 

The "correlation functional," eq. flUf), is recaptured through (FW eqs. 32.14 and 32.16), 

^ = 2Im[^](l-e-^)- 1 . (31) 

Reverting to the single-channel problem for illustrative purposes, we define the polar- 
ization function as 

U(q,uj) = ^[n,n] q ^, (32) 

so that 

S(q, u) = 2Im[n(g, uj)} (1 - e^) -1 . (33) 

Note that we can use the fact that Imil(g, uS) is odd in the variable, lu, to derive from ( |3"3"D 
the relation that embodies detailed balance, S(— uj) = e~^ w S(u). The ring approximation 
(sometimes referred to as the Random Phase Approximation, RPA) then gives 

where v(q) = J d 3 xV(x) and H^°\q, uj) is the free polarization. Eq. (|34]) is of the same general 
form as eq. (0) for the equal-time and small-g limits. If we take only the numerator in 
(IH), we recover the effects of Pauli blocking alone. 

To complete the calculation, we need the real and imaginary parts of H^(q, uj), which 
we derive in Appendix A: 

„ , m/ , m 2 r°°ds, r l + e -( s+( 3) 2 +^ 1 

and 



2tt 

where 



m 2 r 1 4- p-Q +P» n 
ImnW (q,uj) = -V log * R , (36) 



and for m we should substitute the effective mass, m* 130 
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To obtain the various S(q, u) functions that are used in (|34D for determining the rates, the 
above formulae are generalized to the multiple-channel case by considering the correlations 
of every pair of densities from the set n p , n n , n^n^. The cross correlations between the 
densities and the spin densities, nf , vanish, as do all dynamical connections between the two 
at the ring level, so that we are left with two 2x2 problems in solving for the correlations. 
In each sector, we define a matrix polarization function IT(g, u)ij = T\a^ nj\ q>Uj and a ring 
approximation, which is the simple matrix extension of (|3"1|). For the free polarization matrix, 
n(°) we take a diagonal matrix in which the proton and neutron elements are given by (|35|) 
and fl3"rj|), with the respective proton and neutron chemical potentials, /i p and fi n , in place 
of \i. The generalization of the potential, v , will connect the proton and neutron elements. 
The elements of the structure function matrices, S and S A , used in the rate formulae follow 
immediately from the matrix form of fl3"4|). It remains to address the potentials. 

A. Nucleon Potential Model 

A general velocity-independent local two-body interaction between nucleons (a) and (b) 
is given by, 

V a > b = V x {r) + r a -r b V 2 {r) + a a ■ a b V 3 (r) + f a ■ r b a a ■ a b V 4 (r). (38) 

This leads directly to the following construction: [] 
a) to determine the Fermi S(q, u) elements, we use 

Vpp(q) = v nn {q) = vi(q) + v 2 {q) 

v pn {q) = v np (q) = vi(q) - v 2 (q) (39) 
and b) to determine the Gamow-Teller S A (q, uj) elements, we use 

Vp P {q) = v A n (q) = v 3 (q) + v 4 (q) 

= <(?) = V M ~ V M- (40) 

Of course, in the nuclear interaction problem within nuclei, as well as at the higher 
densities in neutron star matter, there are a few complications: 1) there isn't a potential, 2) 
taking one anyway, it is too strong to allow the use of perturbation theory [], 3) the better 
numerical methods (e.g., correlated basis functions) for determining the ground-state energy 



4 In interpretating these matrix elements it should be borne in mind that they are not analogous 
to matrix elements of an operator between single-particle states. They operate in a space of two 
neutral densities (the arguments of J-) and describe the scattering of two particles, with no charge 
exchange. 

5 For example, our plane-wave Hartree approximation, used with the potential that is attributed 
to w-meson exchange, would give much too much positive energy and too great a rate reduction 
for the Fermi terms, since it doesn't keep the particles apart at short distances. 
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of nuclear matter from a phenomenological potential may be poorly adapted to calculation 
of the dynamic structure functions. 

Since it is the excitations of the medium, rather than the cold equation of state, that really 
enter this problem, the Landau Fermi liquid theory provides a framework for proceeding. 
This theory defines an energy functional associated with variations of the various densities. 
We define combinations of density variations, 

80,0 = 8n>p + Sn n 
£3,0 = Sn p - 5n n 
^0,3 = + Snl - 5n l p - Sn^ 

83,3 = £n} - 5nl - 5n l p + 5n l n , (41) 

where the vertical arrows signify spin-up and spin-down densities. The absence of an arrow 
implies the sum of these densities. 

The energy response to these variations that defines the FLT is given by 

5E = Y: ef\q)5n i {q) + ^ BM^M^ + 

q,i q,q> 

<M<j)Mq')*" + M^MqOG + <Mq)<Mq')G"], (42) 

where A = TT 2 {2m*pp)~ 1 and pp is the nucleon Fermi momentum. The FLT parameters are 
usually a function of the angle between q and q'. We shall consider the wavelength in our 
application to be long enough so that only the S-wave parameters, F , Fq, Gq, G' , enter. 
Then, we note that the interaction term, quadratic in the 5s, is exactly what one would get 
from the form (j38|), with zero-range two-body potentials, where we define Vi = J d 3 rVi(r), 
and obtain 



Vl = \F ; v 2 = \Fb v 3 = \G ; v 4 = \G' . (43) 
Note that all the f,s are real. 



The parameter A in (52) implies that our potential has a density dependence. However, 
we view this as simply part of the mechanics of fixing our parameters / d 3 rVi in terms of the 
hodgepodge of nuclear phenomenology, at nuclear densities, and meson-exchange consider- 
ations that went into the fixing of the FLT parameters. We keep the parameter pp fixed at 
its value for nuclear density, and use our derived parameters at all densities. This is clearly 
a sounder procedure at or near nuclear densities, where unextrapolated phenomenology was 
the main input, but we shall calculate the results over a wide region of densities nonetheless. 
We note that the use of Skyrme interactions, expressing an effective energy functional as a 
form quadratic and cubic in the densities, would be subject to the same caveats. 

For supernova and protoneutron star applications, we must include Coulomb forces in 
our interaction Hamiltonian. We can do so by adding to the (zero-range) proton-proton 
force in (^), a Thomas-Fermi screened Coulomb force between protons, 

v pp = J d 3 rV p ' p = v 1 +v 2 ^v 1 +v 2 + 4vre 2 (g 2 + g 2 F )-\ (44) 
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where q TF = 4e 2 7r 1//3 (3fip) 2 / 3 . In the denser regions of the star, the screening momentum is 
larger than the q for typical neutrino scattering. In this case, the Coulomb term in ( pH|) for 
the proton-proton interaction is independent of q and e 2 . f\ 



V. FINAL PROCEDURE FOR CALCULATING v NUCLEON STRUCTURE 

FUNCTIONS 

From eqs. (|17])> (fL9f), (^l|), and (|22"D and taking only the neutron part of the vector-current 
coupling, the differential scattering rate is given by, 



d 2 T 



(4tt 2 )- 1 G 2 h/ E 2 [1 - U(E 2 )} (1 + cos6)(CZyS nn (q,uj 



duj d cos 9 

+ (3 - cos% 2 [S*{q, u) + St(q,u>) - 25£(g,w)], (45) 

where E 2 =Ei — to. 

The structure functions, 5* and S A , are elements of separate 2x2 symmetric matrices. 
For the vector dynamic structure function, S, we have 



\S p Jq,uj) S nn [q,u)J 



' 'pa 

(q,u) S nn (q, 
The structure function matrix is given by, 

S(q, w) = 21m [n^ (q, v) [1 - v(q)U^ (q, uj)}- 1 ] (1 - e^)' 1 (46) 

where 



(llV>(q,u)) 



to 



and nf) and n^ ) are given by the polarization function, defined in (|35|) and (|36|) and eval- 
uated with the proton and neutron chemical potentials, respectively. The potential matrix 



6 The physics of this additional term can best be elucidated by considering for a moment what the 
effects would be if the Fermi term in our weak current coupled only to the protons in the medium 
(rather than almost entirely to the neutrons, as in the standard model.) In this case the scattering 
of neutrinos would be dominated by the scattering from proton-density fluctuations on the order of 
the neutrino wavelength. For long neutrino wavelengths, such fluctuations are strongly suppressed 
by the Coulomb force. When q « qTF the price paid in energy for the fluctuation is measured 
by the increase in the kinetic energies of the neutralizing electrons, and we note that in this limit 
the Coulomb term in (|44| ) is just the second derivative of the electron Fermi energy with respect to 
electron (or proton) density. In the case at hand, where s'm9w ~ 1/2 so that the weak coupling is 
entirely with the neutrons, the Coulomb force between protons turns out nonetheless to be quite 
important, because of the strong coupling of neutrons to protons in the symmetry-energy term. 
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(v x + v 2 + Ane 2 (q 2 + q^ F ) 1 v x -v 2 
V Vi - V 2 V\ + v 2 1 ' 



0;G' 



1.7 



where the v's were defined in terms of Fermi liquid parameters in (|43|). 

In a real calculation, in all the kinematic expressions m is to be replaced by m*. Unfortu- 
nately, the relation of Landau parameters to experimental results depends upon the effective 
mass in model-dependent ways. Taking m* = 0.75m n as our fiducial value for the effective 
mass, we use parameters from references F = — 0.28; Fq = 0.95; Go 

and A = 2.63 x 10" 5 MeV~ 2 , obtaining, 

Vl = -7.4 x lO^MeV -2 
v 2 = 2.5 x 10- 5 MeV~ 2 

^3 = 

v 4 = 4.5 x 10- 5 MeV~ 2 . 



(47) 



For other values of the effective mass, we keep these potentials at the same value, which is 
to say we assume that the Landau parameters are proportional to m*/m. For the dominant 
spin-independent term, v 2 , this accords with the conventional wisdom that the symmetry 
energy per nucleon can be written in the form a((A — Z)/N) 2 . 

The form for the Gamow-Teller matrix, S A (q, u), is the same as that for S, except that 
the potential matrix is replaced by v A 

V A = ( V3 + Vi Vs ~ Va \ 

V V 3 - V 4 V 3 + V 4 / ' 

Taking the matrix inverses leads to the following forms for the combinations of structure 
functions that appear in (f45|) 

S nn (Q, w) = 2Im[n.( 1 °) J D y 1 ](l - e-*")-\ (48) 

(vi + v 2 )U^ - ( Vl - ^) 2 n(°)nf [1 - 47 re 2 (g 2 + q^Tlf - (v, + v 2 )U^]-\ 

(49) 



where 

By = 1 



which corresponds in the (q,ui) —> (0,0) limits to eq. (WA) - 

If, as in (^7|), we take V3 = 0, we obtain the simple result for the axial-current terms, 



21m 



n( )(g, w ) + n(°)(g,^) 



1- Vi [uf\q,u;) + U^(q,cu)] 



(0), 



1-e 



-Pu\-1 



(50) 



which corresponds in the (q, u) — > (0, 0) limits to eq. 



VI. RESULTS 

The formulae that we have developed to calculate dynamic structure functions and scat- 
tering rates (i.e., cross sections) for neutrino-matter scattering in nuclear matter, in partic- 
ular eqs. fl45f) , ([48]), and (|50|) , can now be used to obtain quantitative results. One is free to 
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insert whatever parameters for whatever thermodynamic states and nuclear models, but we 
here choose to focus on a generic subset of possibilities to demonstrate the character of the 
new results. Since there are six neutrino species in thermal equilibrium in the supernova 
core and the electron types have very different chemical potentials than the u^s and u T s, 
we drop the blocking term, (1 — / I/ (£ , 2 )), in calculating the total suppression factors. This 
term is trivial to include in the general case, and for the z/ M s and z/ r s its omission introduces 
only a small error. However, we want to avoid expanding the number of comparisons un- 
duly and the reader is free to employ the derived equations to calculate everything for any 
combination of parameters. Since for the differential cross sections one does not integrate 
over the energy transfer, u, and we present not absolute cross sections, but cross sections 
normalized to the non-interacting case, the differential cross section results are fully general. 
Final-state nucleon blocking is always included. 

For the Fermi term, since Cy = 1/2 — 2sin 2 #iy ~ 0, we drop the proton structure 
function in (IS). Furthermore, we use the potential parameters given in eq. P?|), and in 
eq. ((49|) we drop the third term. This term would have been significant had it not been 
for the Coulomb term in the denominator, an illustration of the importance of the explicit 
inclusion of Coulomb forces, even for the neutron density correlations. Since the i^s are all 
real, we obtain for the structure factors used in (S5j), 

(51) 



S F (q,u) = 2Imnf(l - e-^Cx 



where 



(52) 



and 



where 



S A {q,u) = 2[lmUP(q,u)+lmU^(q, 



to 



c 



A 



1 - v GT (ReIlV\q,Lu) + ReU^(q, u))f + v 2 GT [imllf (q, u) + ImlT^g 



(53) 



(54) 



The F in S F (q,u) and the A in S A (q,u) stand for Fermi and Gamow-Teller (axial) 
and Vp and vqt equal (vi + v 2 ) and t> 4 , respectively. Cy^ is the correction factor due to 
many-body effects for a given momentum transfer (or scattering angle) and energy transfer. 



A. Collective Excitations of the Medium 

For most regions of phase space, C A and are greater than one and represent suppression 
in the scattering rates. Their effects on the integrals over u and 9 are always suppressive. 
However, the terms containing the real parts have roots; these roots represent collective 
excitations. For the Fermi term, zero sound in the medium can be generated if the scattering 
has a (u>,q) pair that satisfies the mode's dispersion relation, i.e., if it hits the resonance. 
Similarly, for the Gamow-Teller term, spin waves in the protons and the neutrons (related 
by a set phase) can be generated. These modes are the traveling-mode equivalents of the 
Gamow-Teller resonance in nuclei (a standing wave). The zero sound of the Fermi part is 
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analogous to the Giant-Dipole resonance in nuclei. The resonances increase the structure 
function when the scattering transfer ratio, u/q, equals the ratio of the collective excitation's 
angular frequency and wavenumber. For a given scattering angle, one can plot the differential 
cross section in u> and cos 8 as a function of u/q to see the resonances. In Figure (0), we 
display this for five different angles between 15° and 180°, an incident neutrino energy of 
20 MeV, a temperature of 5 MeV, a density of 3 x 10 14 gm cm -3 , and an electron fraction, 
Y e , of 0.3. We see in Figure [TJ that the resonances in both the forward and the backward 
directions line up at the same values of ui/q, as expected for a collective mode, and we can 
straightforwardly calculate the mode's dispersion relation. This is akin to the Cerenkov 
effect. Note that the Gamow-Teller term dominates the Fermi term, so that in Figure [I] 
we are really seeing the spin waves related to the Gamow-Teller resonance. However, the 
dispersion relations for zero sound and these spin waves are generally similar. In fact, 
recalling the classic result [35[] that in the weak-coupling limit, the speed of zero sound in 
a degenerate system is ~ Vf ermi , where Vf ermi is the Fermi velocity, and recalling that for 
nucleons in nuclei Vf ermi is ~ 0.3c, the calculated resonance value of u/q is not unexpected. 
In Figure |2|, we plot the Gamow-Teller structure function versus uj/q for various values of 
u, m*, and two values of the density. At m* = m n , for each value of the density we obtain a 
sharp resonance, but at two different speeds, reflecting the crude p^-dependence expected 
for Vf ermi . For a given density, the mode speed is seen in Figure ^| to be inversely proportional 
to the effective mass. The width of the resonance is determined by the magnitude of the 
imaginary part of the polarization function. 



B. Differential Scattering Cross Sections and Suppression Factors 

To calculate the singly-differential scattering cross sections (da /duo) , we must integrate 
([45|) over cos 9. Since for a given uj and incident neutrino energy, E±, this integration is 
also over a range of qs, in the process we are smoothing over resonances. As a consequence, 
there is no obvious direct signature of them in the final result. It is the doubly-differential 
structure functions and cross sections that retain the character of the collective modes and 
resonance. The integral over the singly-differential scattering cross section yields the total 
cross section, and this can be compared to that without correlations to gauge the magnitude 
of the suppression of the rate. To demonstrate the nature of the correlation effects, we have 
opted to present figures and tables for a subset of the possible p-T-Y e -E\ combinations. 
Recall that for this purpose we have employed the default potentials (^) and effective mass 
(m* = 0.75m n ). Note, however, that the formalism we have derived is fully general, as 
long as the nucleons can be assumed to be non-relativistic and we are in the perturbative 
limit. Since we calculate at given neutron and proton densities, p n and p p are implicit in 



the formalism, as would be the shifts in them due to interactions |30 



We depict in Figure [3] the singly-differential scattering cross section, divided (normalized) 
by the total scattering cross section off of nucleons (a\), for a range of incident neutrino 
energies. We have used eqs. (|45|), (|5T|), (|52"D, (|53"D, and fl54|) to generate these curves. Figure |3| 



is a study of the dependence on incident neutrino energy (1, 5, 10, 20, 30 MeV), at a fiducial 
(but arbitrary) density (p = 3 x 10 14 gm cm -3 ), temperature (T= 5 MeV), and electron 
fraction (Y e = 0.3). Also shown as dashed lines are the corresponding curves without the 
many-body correlation effects (Cy^) and for m* = m n , but with final-state nucleon blocking. 
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Table I depicts the corresponding many-body suppression factors, obtained by integrating 
under the curves, as well as both the average energy transfers (<u >) and the rms of the 
energy transfers (uj rms ). Without neutrino blocking, the total suppressions we present in 
the tables are very close to the full results for u^s and u T s. For v e s, it is more important to 
include the (1 — f^E?)) term, since there is a net electron lepton number in protoneutron 
stars for most of their interesting lives. Figure |3| is rich with information that we will try to 
summarize. Positive uj represents energy lost from the neutrino and negative u represents 
energy gained by the medium. First, notice that the widths of the curves increase with 
the incident neutrino energy. Even for the curves without many-body effects, the widths 
as depicted in Figure and in Table I are quite large. Such widths call into question the 



elastic approximation [|26| , |28| , but also imply that energy equilibration for the v^s and v T s 
by z/-nucleon scattering dominates over the annihilation process, v + v — > e + + e~, and 
^-electron scattering, both with much lower cross sections | 29fl . Equilibration for the z/ e s is 
still via the charged-current absorption process, v e + n — > e~ + p. 

One major reason the widths are larger than are familiar is that in the past people 
thought that the neutrino could lose in i/-nucleon scattering an energy equal to only about 
—Ef/m n c 2 , i.e. that the fractional energy lost is of order p u /m n c (~1%). However, this 
assumes that the nucleons are stationary. In fact, they are thermal and, the fractional energy 
they can transfer in a collision to the neutrino is of order p n /m n c. Since the nucleons have 
such a large mass, if they and the neutrino have the same energy, p n /m n c is much larger than 
p u /m n c, at incident neutrino energies of 10-30 MeV by as much as an order of magnitude. 
The formalism we employ incorporates the kinematics of such a collision, a realistic Fermi- 
Dirac energy distribution for the nucleons, and final-state nucleon blocking. The upshot is 
the broad distributions, even without the Cy^s , seen in Figure |3| and tabulated in Table 
I. Including many-body effects further flattens and broadens the distribution (see below), 
while lowering the central values of da/du, as well as the total integral over u. 

As Table I demonstrates, even without many-body effects, nucleon blocking is a large 
effect, larger than the elastic |— correction that comes from the low-temperature expansion 
in powers of T. That correction can be applied in the extreme degenerate limit, but the 
nucleons are only partially degenerate in protoneutron stars and supernovae, even at nuclear 
densities, and the extreme degenerate limit is never achieved. The error in using the |— - 
correction ansatz can be as much as a factor of two, depending upon the nuclear interaction 
model. 

Note in Table I that the many-body cross section suppression factors are not strong 
functions of the incident neutrino energy, except at very low energies (e.g., 1 - 5 MeV), 
but that the suppression factors themselves are quite large at this fiducial thermodynamic 
point (near nuclear density). The extra suppression, beyond that due to nucleon blocking, 
is due both to the decrease in the effective mass (factor of ~1.5) and to the many-body 
correlation terms, Cy^ • The total factor is between 10 and 30. This means that the 
combination of final-state nucleon blocking (a factor of ~6 by itself) and many-body effects 
(another factor of 3 to 4) renders the supernova core much more transparent to neutrinos 
than previously thought. We remind the reader that we are here calculating only the neutral- 
current rates. While they dominate for the z/ M s and u T s, the charged-current dominates for 
the v e s. Nevertheless, we are in the midst of preliminary calculations that indicate that the 
suppression of these rates is also quite large (Burrows & Sawyer, in preparation; see also 
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Reddy et al, in preparation). 

The previous assumptions concerning the distribution in energy transfers for z/-nucleon 
scattering were more akin to the narrow (quasi-delta function) curve, seen in Figure |3] for 
1-MeV incident neutrino energy. The contrast between that curve and the others is manifest. 

Figure f| is similar to Figure [| but for a range of temperatures (5, 7, 10, 15, 20, 30 
MeV). The calculations were done at E\ = 20MeV. Superposed are the corresponding 
reference curves for C^,a — 1 an d rn* = m n . Table II lists the suppression factors, < uj >, 
and uj rms . As the temperature rises, the suppression diminishes. It is clear in Figure |] that 
energy transfer from the medium to the neutrino becomes more likely as the temperature 
rises. This is to be expected and is all the more pronounced in the many-body case. As 
in Figure ||, the differential cross section is flattened when interaction effects are included, 
but Figure f| shows that this flattening effect is more significant at higher temperatures. As 
Table II demonstrates, though the total suppression diminishes with temperature, many- 
body effects still increase it, the more so at lower temperatures. Given the large values of 
< uj > and uj rms > h is difficult to see how even a Fokker-Planck treatment of neutrino energy 



redistribution via these processes could be viable and the full redistribution formalism PD|J^1 
may be required, at least for the v^s and the v T s. 

Figure |5| depicts the density dependence of the singly-differential cross section from 
p = 10 12 gm cm~ 3 to p = 3 x 10 14 gm cm -3 , with and without many-body effects. The 
temperature is kept constant at 5 MeV and the incident neutrino energy is 20 MeV. Table 
III displays the corresponding total suppression factors, < uj >, and uj rrns . Both Figure || 
and Table III demonstrate that the total suppression effect increases quickly with density 
beyond 10 13 gm cm -3 and that that due solely to many-body effects increases quickly 
beyond 10 14 gm cm -3 . In fact, the total correlation suppression factor reaches ~100 at 10 15 
gm cm -3 , calling into question the perturbative assumption itself. Nevertheless, it is clear 
that correlation and many-body effects can drastically lower the neutral-current rates in 
high-density supernova cores. Note that the widths displayed in Figure [5] increase with 
increasing p. This reflects the increasing nucleon degeneracy and increasing average nucleon 
energy that accompanys increasing p. 

Taking a temperature, density and Y e profile from an early post-bounce model of Bur- 
rows, Hayes, & Fryxell ||, we show in Figure |] the corresponding differential cross section 
curves at E\ = 20 MeV for a "realistic" profile. Table IV displays the total suppression 
factors, <uj>, and uj rms . The total suppression effect at the center is greater than 30, that 
due solely to many-body effects is ~5. The <uj> and uj rms are also large. Even at 10 14 gm 
cm~ 3 , the extra suppression effect due to many-body effects is ~2. 

It is interesting to compare the suppression factors calculated from our full equations 
with those that come from the long-wavelength limits given in eqs. ([27]) and (|29|), using the 
same potential parameters. At a density of 3.0 x 10 14 gm cm -3 , a temperature of 5 MeV, 
and for neutrino energies less than 10 MeV, the factors, calculated the two different ways, 
are within 10% of each other. However, for a neutrino energy of 30 MeV, the full calculation 
gives a suppression that is 30% greater than in the long- wavelength limit. It will be recalled 
that our "long-wavelength" limit also requires that the energy transfer, uj, be small. Since 
the energy transfers are of the order of PfQ, this limit can be achieved with fixed nucleon 
mass only by going to lower density. With temperature fixed at 5 MeV, and neutrino energy 
fixed at 20 MeV, we find less than a 5% discrepancy for a density of 10 13 gm cm -3 . However, 
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the discrepancy rises steadily to nearly 50% at a density of 10 14 gm cm -3 . These numbers 
confirm that the static approach is adequate only in some of parameter realms occupied by 
super novae. 



C. Correction Factors versus Energy Transfer 

It is instructive to calculate d 2 T / duj d cos 9 , after integration over cos#, {i.e., da/duj), 
with and without the Cy and C4 terms and a renormalized nucleon mass, and to take the 
ratio. This gives one a sense of the integrated correction factor to the differental rate due to 
many-body effects. This is the many-body correction factor to what might be called S(ui). 
Figures [7|, H and |I] depict the logarithms of these correction factors as a function of to. 
These figures correspond to the Figures | through ||. We see the expected suppression factor 
in the small-u; regime that dominates the total suppression integral, but at high \u\s we see 
manifestations of the resonances. The rapid decrease in ImlT )(g, lo) with |u;| defeats the 
increase in the correction factor at high \u\ in the total cross section integral. However, this 
resonance effect is partly responsible for the high— \u\ flattening we see in Figures^ through 
|6|. Note that this correction factor does not go to zero at u = 0, can go far above 1 for high 
\w\s, and does not have the form, u 2 /(u 2 + K). 



VII. AN ILLUSTRATIVE MODEL OF PROTONEUTRON STAR COOLING 

While the consequences of the inclusion of many-body effects on neutrino cross sections 
and energy redistribution are by no means clear, one can straightforewardly ascertain the 
potential of such effects to alter aspects of supernova and protoneutron star development 
by performing simplified protoneutron star cooling calculations. We have chosen to inves- 
tigative the evolution of the electron neutrino luminosity (L Ue ), under certain simplifying 
assumptions, with and without various ad hoc alterations to the neutrino opacities. The 
code, developed by Burrows & Lattimer ||, uses the Henyey technique with a simple nu- 
clear EOS, is general-relativistic, and handles the transport of neutrinos of all species in the 
diffusion approximation. It conserves total energy to about one percent over 20 seconds of 
evolution. The default cases are models of L Ue versus time after bounce without accretion 
and with an assumed accretion rate of O.4M e"^ 7 "", where r a is 0.5 seconds. The initial 
baryon mass is 1.3 M Q . We alter the total opacities in an ad hoc fashion to mimic the de- 
creases we can anticipate given our preliminary exploration of neutral-current suppressions 
and compare the results to the default models. The two exploratory models we present 
assume that the total opacities are decreased by a fixed amount above a given density, in 
one case 5 x 10 13 gm cm -3 and in the other 10 14 gm cm -3 . For the 5 x 10 13 gm cm -3 
case the suppressions were 0.3, 0.1, and 0.05 (extreme) and for the 10 14 gm cm -3 case the 
suppressions were 0.3 and 0.1. These were guided by our results, but should be considered 
arbitrary. Note that the more conservative case for which the opacity is altered modestly 
and at the higher density is the more likely. 



The results are plotted in Figure |TT|. From a comparison between the fiducial model 
and the one without accretion, one notes that accretion dominates as a power source in the 
early seconds. If the supernova is reignited within the first hundreds of milliseconds, it is 
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unlikely that the new opacities at high densities will play a central role in its revitalization. 
However, if the delay to explosion is many hundreds of milliseconds to seconds, higher L Ve s 
occasioned by the many-body and final-state nucleon blocking effects may well play a role in 
the supernova, and in powering the explosion after it is relaunched. Even partially enhanced 
luminosities have been shown to be important for neutrino-driven explosions . As 

Figure [TT] demonstrates, the lower the density at which opacities are altered, the earlier 
the effect. The enhancements in the luminosities after a few seconds can be a factor of 
two, though enhancements of 30%-50% are more likely. Since energy is conserved, the 
luminosities at later times (5-60 seconds, not shown) are decreased relative to the those of 
the fiducial model. Within this model set, the theoretical signals in Kamioka II || and 1MB 
range by no more than a factor of two. Drastic decreases in the total opacities translate 
into only moderate (but intriguing) changes in L Ue and the signals because even a completely 
transparent core merely flattens the interior temperature gradients. The energy must still 
diffuse through the unaltered density region below 5 x 10 13 or 10 14 gm cm -3 . 

We point out that for progenitor stars in the 8 M to 13 M mass range early accretion 
onto the protoneutron star before explosion will be slight. Hence, alterations in the high- 
density neutrino cross sections might show themselves earlier in lower-mass progenitors, in 
time to be more unambiguously of importance in reigniting and powering the supernova. 
Nevertheless, the reader should be cautioned that these conclusions await a more thought- 
ful implementation of these new neutrino opacities and that the calculations depicted in 
Figure O are merely suggestive. 



VIII. CONCLUSIONS 

In this paper, we have developed a consistent formalism for calculating the Fermi and 
Gamow-Teller structure functions in nuclear matter, including many-body correlations and 
the full kinematics of neutrino-nucleon scattering. Previously all but ignored, the effects 
of density-density and spin-spin correlations on neutrino scattering rates are found to be 
significant. Above nuclear density, the total cross section suppression factor due to final- 
state nucleon blocking, correlations, and a reduction in the nucleon effective mass can be 
more than an order of magnitude. The upshot is that supernova cores are more transparent 
than previously thought. Since the many-body corrections to the charged-current rates for 
the z/ e s have yet to be published, in this paper we conclude only that the total opacities 
of nuclear matter to z/ M s and z/ r s (and their anti-particles) are qualitatively altered. The 
magnitude of the rate corrections due to many-body effects increases with density, decreases 
with temperature, and is roughly independent of incident neutrino energy. 

In addition, we find that the neutrino-matter energy transfer rates due to neutrino- 
nucleon scattering are much larger than heretofore assumed. As a consequence, energy 
equilibration for the non-electron-type neutrino species may well be by what was previously 
considered an elastic process. This obtains even for free nucleons at high densities, but is all 
the more true with many-body effects included. An identifiable part of the energy transfer 
enhancement comes from the excitation of collective modes in the matter, modes akin to the 
Gamow-Teller and Giant Dipole resonances in nuclei. Energy is transferred by the Cerenkov 
mechanism to or from these modes when the energy transfer (±u;) and momentum transfer 
(q) satisfy the dispersion relation of the medium's excitations (zero sound and spin sound, 
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both of which require repulsive interactions). 

Under the assumption that both the charged-current and the neutral- current cross 
sections are decreased, we calculated a set of ad hoc protoneutron star cooling models 
to gauge the importance of the new opacities to the supernova itself. While the early 
luminosities are not altered, the luminosities after many hundreds of milliseconds to seconds 
can be altered by factors that range from 10% to 100%. Such factors may have a bearing on 
the efficacy of the neutrino-driven supernova mechanism, the delay to explosion, the energy 
of the explosion, and the strength and relative role of convective overturn. The magnitude 
of the relative enhancement of the driving luminosity is a function of the post-bounce mass 
accretion rates and may be larger for the least massive, massive stars. Hence, these new 
opacities may be more germane to the terminal behavior of the massive stars most favored 
by the IMF. 

It remains to derive the effects of correlations on the charged-current rates, as well as to 
develop a simple algorithm for incorporating these effects into supernova and protoneutron 
star codes, with full energy redistribution. Since many-body effects at the highest densities 
are quite large, a non-perturbative approach may need to be developed for the deep interiors. 
Furthermore, the structure function calculations should be done in a manner fully consistent 
with the equation of state employed, since the physics of the two are inextricably linked. 

If we find that there is indeed an across-the-board reduction in the neutrino-matter 
rates at high densities, it will be yet another reminder that the keys to the supernova puzzle 
and its systematics lie not in any one realm of expertise (radiative transfer, hydrodynamics, 
the equation of state, the weak interaction, etc.), but in all, and that the 20% effects we 
seek under lamp posts might at times divert us from discovering those far larger just a few 
steps down the street. 



ACKNOWLEDGMENTS 

We would like to acknowledge stimulating and productive conversations with S. Reddy, 
M. Prakash, G. Raffelt, D. Seckel, and J. Lattimer, as well as support from the NSF under 
grant No. AST-96-17494. We would also like to acknowledge the germinating influence of 
the Santa Barbara Institute for Theoretical Physics, supported by the NSF under grant No. 
PHY94-07194. 



APPENDIX A 

To get expressions for the real and imaginary parts of that are efficient for compu- 
tation, we begin with FW (33.4), with notations changed to those of this paper 

rf°\q,u) = -2 / /p+q ~ /p . (55) 

1 ' J (27r) 3 ^ + ^ + e p -e p+q 1 ; 

We can write each of the Fermi-Dirac distributions that appear here as a sum of Boltz- 
mann distributions, 

oo 

/(p) = £(-iy + V^). (56) 

3=1 
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We use this expansion, together with FW (33.8) for the imaginary part of IT ) in the 
Boltzmann limit, to obtain, 

m 2 oo 

Imn(%o;) = ^J2(-iy+^[e^-Q 2 ) - e i(^-Q 2 -^)] ) (5?) 

where 

Q=( ^ ( _- + X) (58 ) 
2 q 2m 

and the terms with j > 1 in the series come from the replacement, (3 —>■ j/3, in the Boltzmann 
result. The sum is easily performed to give the result shown in fl3~6"|). For the real part, we 
use the same trick, beginning with the Boltzmann form (FW 33.9), 

2 

ReU^(q,cu) boltz = -^— e -P^{Q) + (u - -u), (59) 

where 



$(Q) =2Q [ l dye Q ^ 2 - l \ (60) 
Jo 

We obtain, 

R*n (0) (?>") = T^rQ / dyEC- 1 )^'- 1 ^^^-^ + (w - -u,). (61) 
2i\ 6 i l \3q Jo pr{ 

Next, we represent the j -1 ^ 2 factor under the sum by 7T- 1 / 2 JZo ds e-i s \ so that we can 
sum the geometric series. Displacing the s integration variable s — > s + QY, we obtain, 

m 2 ^oo /•! e -(s 2 +2sQj/+Q 2 +/3 At ) 



The ?/ integral is now easily performed to give the result (|35f) . 



APPENDIX B 

We show directly that in the large-m, small-g limits the complete correlation functions as 
calculated in the ring sum in §IV reduce to the long-wavelength limits that were determined 
from the single-particle energies in §11. In the large-m limit, the function Imn°(g, uj) becomes 
more and more concentrated around uj = 0, so that to calculate the rate we set uj equal to 
zero, except in the factor, 

/oo 
dujlmli(q,uj)(l- e'^y 1 . (63) 
-oo 

In the limit, we make the replacement (1 — e - ^) — > (3uj and use the Kramers-Kronig 
dispersion relation in the variable uj to write I\ — > ReII(g, uj = 0). Finally, we use fl3"oD in 
the limit Q — > to show that ReII^(0,0) = —h(/i), where h(fx) is the function defined in 
(0). The ring approximation for the one-channel case ( |54"D with a potential U now reads, 

h^h(ji)/(l + Uh(fi)), (64) 



in agreement with (|14j) 
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TABLES 



TABLE I. The total suppression factors (S), average energy transfers (<cj>) and rms energy 
transfers (w rms ) for various incident neutrino energies (E\). The temperature is 5 MeV, the den- 
sity is 3 x 10 14 gm cm~ 3 , and the electron fraction is 0.3. These quantities are shown for both 
the many-body case and for the case without correlations or a renormalized mass (subscript or 
superscript 0). All the energies are in MeV. For the many-body case, the default nuclear model 
described in the text is employed. The suppression factors are the factors by which the default 
total cross section should be multiplied to obtain the corrected total cross section. 



E 1 (MeV) 


S 


So 


<lo> (MeV) 


1 


0.131 


0.217 


-0.917 


5 


0.090 


0.202 


-3.749 


10 


0.060 


0.186 


-5.064 


15 


0.045 


0.173 


-4.204 


20 


0.038 


0.164 


-2.156 


25 


0.035 


0.158 


0.287 


30 


0.034 


0.154 


2.779 


40 


0.035 


0.152 


7.583 


50 


0.038 


0.155 


12.109 


60 


0.041 


0.161 


16.413 


70 


0.046 


0.169 


20.549 



<uj> (MeV) uo rms (MeV) co° rms (MeV) 



-0.207 0.592 0.295 

-0.852 3.088 1.432 

-3.000 6.068 2.755 

-1.288 8.241 3.939 

-0.970 9.492 4.977 

-0.385 10.195 5.881 

0.402 10.687 6.671 

2.375 11.569 8.007 

4.649 12.522 9.143 

7.063 13.574 10.177 

9.532 14.712 11.161 



TABLE II. As in Table I, the total suppression factors (5), average energy transfers (< uo >) and 
rms energy transfers (ay ms ), but for various temperatures at a given incident neutrino energy (E\ 
= 20 MeV). The subscript or superscript corresponds to the case without many-body effects, but 
with final-state nucleon blocking. This table corresponds to Figure 4, but also includes numbers 
for T = 3 MeV. 



T (MeV) 


S 


So 


<co> (MeV) 


<co>o (MeV) 


L0 rrns (MeV) 


uj° rms (MeV) 


3 


0.021 


0.093 


2.693 


0.598 


6.623 


4.292 


5 


0.038 


0.164 


-2.158 


-0.970 


9.497 


4.977 


7 


0.065 


0.240 


-6.153 


-1.856 


11.221 


5.358 


10 


0.118 


0.350 


-10.193 


-2.699 


12.482 


5.810 


15 


0.227 


0.516 


-14.355 


-3.734 


13.928 


6.614 


20 


0.351 


0.659 


-17.572 


-4.726 


15.695 


7.618 


30 


0.639 


0.900 


-24.603 


-7.118 


21.660 


10.575 
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TABLE III. As in Tables I & II, the total suppression factors (5), average energy transfers 
(< lo >) and rms energy transfers ( w rms ), but for various mass densities (p) at a given incident 
neutrino energy {E\ = 20 MeV) and for a temperature of 5 MeV. The subscript or superscript 
corresponds to the case without many-body effects, but with final-state nucleon blocking. This 
table is related to Figure 5. 



p (gm cm 3 ) 


S 


So 


<uj> (MeV) 


10 12 


0.897 


0.952 


-0.295 


10 13 


0.560 


0.794 


-0.455 


3 x 10 13 


0.307 


0.591 


-0.766 


10 14 


0.117 


0.162 


-1.525 


3 x 10 14 


0.038 


0.164 


-2.156 



<o;> (MeV) uj rms (MeY) w r ° ms (MeV) 



0.216 2.556 2.157 

0.256 3.211 2.359 

-0.337 4.269 2.725 

-0.562 6.532 3.598 

-0.970 9.492 4.977 



TABLE IV. As in the other tables, the total suppression factors (5), average energy transfers 
(< to >) and rms energy transfers (w rms ), but for a post-bounce profile from Burrows, Hayes, & 
Fryxell (1995), at a given incident neutrino energy {E\ = 20 MeV). The subscript or superscript 
corresponds to the case without many-body effects, but with final-state nucleon blocking. This 
table is related to Figure 6. The density and temperature at the edited points in the model are 
indicated. 



p (gm cm 3 ) 


T (MeV) 


S 


So 


<lo> (MeV) 


<uj> (MeV) 


LV rms (MeV) 


uj° rms (Me 


3.945 x 10 14 


5.053 


0.028 


0.136 


-2.142 


-1.136 


10.082 


5.426 


3.084 x 10 14 


5.518 


0.043 


0.180 


-3.305 


-1.261 


10.135 


5.137 


2.049 x 10 14 


9.457 


0.147 


0.402 


-7.322 


-2.108 


10.290 


5.130 


1.281 x 10 14 


10.561 


0.243 


0.547 


-5.797 


-1.882 


8.656 


4.706 


6.463 x 10 13 


13.531 


0.477 


0.792 


-4.629 


-1.946 


7.460 


4.659 


4.322 x 10 13 


14.554 


0.608 


0.880 


-4.065 


-1.959 


6.994 


4.666 


2.669 x 10 13 


15.611 


0.750 


0.953 


-3.660 


-2.013 


6.676 


4.735 


5.972 x 10 12 


11.287 


0.876 


0.968 


-1.740 


-1.159 


4.554 


3.609 


1.082 x 10 12 


6.137 


0.907 


0.957 


-0.509 


-0.369 


2.874 


2.416 
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FIGURES 




-.6 -.4 -.2 .2 .4 .6 



w/q 

FIG. 1. Logio °f the doubly-differential cross section for neutral-current neutrino-nucleon 
scattering versus u>/q for scattering angles 15°, 45°, 90°, 135°, and 180° . The calculations were 
performed at a temperature of 5 MeV, a Y e of 0.3, a p of 3 x 10 14 gm cm -3 , and an incident 
neutrino energy of 20 MeV. The default potentials (vqt = 4.5 x 10~ 5 and vf = 1-76 x 10 -5 ) and 
effective mass (m* = 0.75 m n ) were employed. The differential cross section is divided by the total 
scattering cross section (<7i) in the non-interacting, no-nucleon-blocking, u = limit. 
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8 



p = 3x10 , 10 gm cm 

T = 5 MeV; Y e = 0.3 

e„ = 20 MeV; u = 6, 8, 10 MeV 

v GT = 4.5xl0 -5 ; m* = 1.0, 0.75 m n 




.2 .3 .4 .5 .6 .7 .8 

w/q 



FIG. 2. Logio °f the Gamow-Teller structure function versus ui/q for an incident neutrino 
energy of 20 MeV, energy transfers, ui, of 6, 8, and 10 MeV, two values of the effective mass 
(m* = [0.75m n , 1.0m n ]) and two values of the density (p = 3 x 10 14 and 10 15 gm cm" 3 ). A 
temperature of 5 MeV and a Y e of 0.3 were used, as was the default vqt (= 4.5 x 10~ 5 ). 
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(j (in MeV) 

FIG. 3. The logio of the differential cross section for z/-nucleon scattering versus the energy 
transfer, co, for various values of the incident neutrino energy {e v = 1, 5, 10, 20, 30 MeV). The 
dashed curves neglect the many-body effects associated with m* and Cy,.A) while the solid curves 
include them. A density of 3 x 10 14 gm cm -3 , a temperature of 5 MeV, and an electron fraction, 
Y e , of 0.3 were assumed. The curves were normalized to the total i/-nucleon scattering cross section 
without nucleon blocking or many-body effects. 
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a) (in MeV) 



FIG. 4. Similar to Figure 3, but for various temperatures (5, 7, 10, 15, 20, 30 MeV) and at an 
incident neutrino energy of 20 MeV. 
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T 




a) (in MeV) 



FIG. 5. Similar to Figure 3, but for various densities (10 12 , 10 13 , 3 x 10 13 , 10 14 , and 3 x 10 
gm cm' 3 ), a temperature of 5 MeV, and an incident neutrino energy of 20 MeV. 



33 



i i i r 




a) (in MeV) 



FIG. 6. Similar to Figure 3, but for an early post-bounce temperature-density- Y e profile from 
Burrows, Hayes, &; Fryxell (1995) (see Table IV). The incident neutrino energy was assumed to be 
20 MeV. 
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FIG. 7. Logio of the correction factor due to both Cv,.4 and an effective mass of 0.75m n that 
corresponds to Figure 3 (the e„ = 1 MeV line has been omitted). There is suppression at low \uj\s, 
but a resonant enhancement at high |o>|s. Note that at very small |w|s the factor is a weak function 
of incident neutrino energy. 
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FIG. 8. Logio of the correction factor due to both Cv,a an d an effective mass of 0.75m n that 
corresponds to Figure 4. There is suppression at low \co\s, but a resonant enhancement at high \ui\s. 
Note that at small |u;|s the magnitude of the correction is a decreasing function of temperature. 
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w (in MeV) 



FIG. 9. Logio of the correction factor due to both Cv,.4 and an effective mass of 0.75m n that 
corresponds to Figure 5. There is suppression at low \uj\s, but a resonant enhancement at high 
\oj\s. Note that at small \uo\s the factor is a strong function of density. 
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FIG. 10. Logio of the correction factor due to both C\j,a an d an effective mass of 0.75m n that 
corresponds to Figure 6. There is suppression at low |u;|s, but a resonant enhancement at high 
\uj\s. Note that at small |w|s the magnitude of the correction is a strong function of position in the 
star and is largest at the center. 
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500 1000 1500 2000 2500 

t (in milliseconds) 

FIG. 11. Logio of the electron neutrino luminosity (L Ue ) in units of 10 51 ergs s _1 versus time 
after bounce in milliseconds, with and without accretion. For the accretion models, total opacity 
suppression factors of 0.3, 0.1, and 0.05 were assumed above 5 x 10 13 gm cm" 3 and of 0.3 and 0.1 
were assumed above 10 14 gm cm -3 . The fiducial model is dashed, the model without accretion 
is dot-dashed, the models with correction above 5 x 10 13 gm cm -3 are dotted, and those with 
correction above 10 14 gm cm -3 are solid. On this plot, the models with the largest corrections 
have the highest luminosities after 2500 milliseconds. The comparisons between the dashed curve 
and all others are the most germane. 



39 



